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Abstract. We solve the boson normal ordering problem for (9(0^)0 + w(at))" with 
arbitrary functions q and v and integer n, where a and are boson annihilation 
and creation operators, satisfying [a,a^] — 1. This leads to exponential operators 
generalizing the shift operator and we show that their action can be expressed in 
terms of substitutions. Our solution is naturally related through the coherent state 
representation to the exponential generating functions of Sheffer-type polynomials. 
This in turn opens a vast arena of combinatorial methodology which is applied to 
boson normal ordering and illustrated by a few examples. 



1. Introduction 

Three decades ago Navon [1] and Katriel [2] observed that the normal ordering of 
quantum operators is a combinatorial problem. Navon considered general products 
of fermion creation and annihilation a operators and showed that such monomials 
can be expressed as sums of lower order normal products with coefficients which are 
combinatorial numbers called rook numbers. Katriel analysed the normally ordered 
form of (a^a)"' for the boson case a,a^ =1 and arrived at a solution given in terms of 
the Stirling numbers of the second kind. Following these ideas we presented closed form 
solutions for the normally ordered form of Taylor expandable boson operator functions 
F{w{a, a"!")) where w{a, a^) is either a monomial in a and operators or a monomial in 
a-|-a^ [3]. Using methods of combinatorial analysis we identified combinatorial numbers 
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appearing in such normally ordered expressions and described their properties. Recently 
Witschel [4] presented an alternative technique for this problem. 

In this Letter we solve the problem of normally ordering expressions of the form 
{q{a'')a + v{a''))"' for arbitrary functions ^(a^) and v{a''). This generalizes results [5], [6] 
obtained in the mid 80 's. Our considerations have been inspired by operational methods 
proposed by Dattoli et al. [7], [8]. Our approach connects the boson normal ordering 
problem to the notion of substitution groups and, if one employs the coherent state 
representation, to the theory of Sheffer-type polynomials appearing in various areas of 
combinatorics [9], [10]. 



2. Operational formulas via substitution group 



We consider bosons satisfying a, = 1 , and use the representation of the Heiscnbcrg 
- Weyl algebra in terms of multiplication and derivative operators 
with the Taylor formula exp (A-r-j F{x) = F(x + A) and consider a generalization of 



-f.x 

ax ' 



1. We start 



the shift operator of the form Egy{X) = exp A {q{x)-^ + v(x)j , where q{x) and v{x) 
are arbitrary functions. We shall find a formula for Eg^y{X)F{x). In fact, as we will 
demonstrate, the following equality holds 



exp 



A {q{x)£ + v{x))]f{x) = g{X, x) ■ F{T{X, x)) 



where 



^^k^^q(T(X,x)), 



dg{X,x) 

dx 



v{T{X,x))-g{X,x) 



T{0,x) = X , 
g{0,x)^l. 



(1) 

(2) 
(3) 



We now exploit basic properties of Eqs.(l)-(3). First observe from Eq.(l) that the action 
of Eq^y{X) on a function F{x) amounts to: a) a change of argument x — > ^(A, x) in F{x) 
which is in fact a substitution; b) multiplication by a prefactor ^(A, x) which we call a 
prefunction. We also see from Eq.(3) that g{X,x) — 1 for v{x) = 0. Finally, note that 
Eq^vW with A real generates an abelian, one-parameter group, implemented by Eq.(l); 
this gives the following group composition law for ^(A, x) and g{X, x): 

T{X + e,x))=T{e.T{X.x)), 

g{X + e,x)^g{\x)-g{e,T{X,x)). ^ ' 

In order to prove Eqs.(l)-(3) first recall the exponential mapping formula which for 
operators A and B states that 



gAA^g-AA ^ 5 ^ y^^j^^ ^] ^ [A, B]] + 



(5) 
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Calculating appropriate commutators one can check that for arbitrary functions v{x) 
and G{x), where x is the multiplication operator, the following operator relation holdsf 

e4£+"W]G'(x)e-4^fe+-(-)] = G(x + A). (6) 

Next, in the l.h.s. of Eq.(l) let us change the variables according to a; = f{y) subject 
to a constraint 

if{y) = q{hy)), (7) 

which leads to 

We treat the r.h.s of above equation as an operator acting on a constant function 1 and 
we get 

which using Eq.(6) gives 

eM£+'^W]F(a;) = F(/(y + A))^(A, y) (10) 

where 

5(A,y) = e4i^+-(/M 1 . (11) 
We can find the differential equation satisfied by ^(A, y) 



-(/(3/))]^(/(^))e-A[^+.(/(2/))]gA[^+.(/(,))] ^ (12) 



= vifiy) + XrgiX,y), 

with the use of Eq.(6) and definition Eq.(ll). The initial condition is g{0,y) = 1. 
If we identify 

T{X,x)=f{f-\x) + X) 

giX,x) = ~giXj-\x)) ^ ' 

then Eqs.(7), (10) and (12) give Eqs.(l)-(3) and so complete the proof. Note that 
formulas of the type Eq.(13) were used by G.A. Goldin in his investigations of current 
algebras [11] and more recently by G. Dattoli et al [8]. For a geometrical interpretation 
of Eq.(13) and the existence of the prefunction see also [14]. 

Here we list several apphcations of Eqs.(l)-(3) for some choices of q{x) and v{x). 
Since Eqs.(2) and (3) are first order linear differential equations we shall simply write 

t A word of explanation is appropriate when comparing Eqs.(l) and (6): Eq.(l) is an equality between 
functions, whereas Eq.(6) is an operator equation. Upon rewriting Eq.(6) as e^[^+''(^)]G(a;) = 
G{x + A)e^[^+^("')] and by acting with it on a constant function 1 we obtain a particular case of 
Eq.(l) for q{x) = 1 and g{X,x) = e^[^"'""^'^^] 1. For the equation for g{X,x) see below. 
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down their solutions without dweUing on details. First we treat the case of v{x) — 0, 
which implies g{X,x) = 1: 

Ex.1: qix)^x, T{X,x)^xe^ (14) 

which gives exp (Ax^j F{x) = F[xe^), a well known illustration of the Euler dilation 



operator exp {Xx^^ . 

Ex.2 : 



q{x) — x^, r > 1, T{X,x) 



X 



(l-X{r-l)x'-')-' 

The above examples were already considered in the literature [3], [7], [12], [14]. 
We shall go on to examples of v{x) leading to nontrivial prefunctions: 

Ex.3 : q{x) = 1; v{x) — arbitrary, 

r(A, x) = X + X, 



g{X, x) = exp 



Ex.4 : 



J du v{x + u) 



2. 



(16) 



Ex.5 : 



q{x) = x; v{x) — x ; 
T(A, x) = xe^, 
g{X, x) = exp Y (^^^ ~ 1 

q{x) — x^, r > 1, v{x) — x^; 
T(X,x)- 



[17) 



(l - A (r - 1) x""^) 



g{X,x) = exp 



x 



s—r+l 



T 



^(l-A(r-l)x"-')''^-^' 



- 1 



(18) 



Closer look at above examples (or at any other example which the reader could easily 
construct) indicates that the group property Eq.(4) may be true only if A satisfies 
certain restrictions arising from the requirement that ^(A, x) and ^(A, x) remain real- 
valued functions. Evidently any result obtained with Eq.(l)-(3) should be examined in 
this respect and resulting restrictions on A be kept in mind. In general we can say that 
the group property Eq.(4) may be vahd only locally [13], [14]. 



3. Implications for boson normal ordering 

The results elaborated above will be used now to treat the problem of the normal 
ordering of operator functions of boson operators. For a general function F{a, a^) its 
normally ordered form J\f ^F(a, a^) = F(a, a^) is obtained by moving all the annihilation 
operators a to the right, using the commutation relations. We may additionally define 
the operation : G{a,a^) : which means normally order G{a,a^) without taking into 
account the commutation relations. Using the latter operation the normal ordering 
problem is solved for F{a, a^) if we are able to find an operator G{a, a' ) for which 
F(a, a^) =: G{a, a^) : is satisfied. To obtain the normally ordered form of functions of 
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boson operators is a nontrivial task even for simple functions. But it is nevertheless 
worthwhile, because the knowledge of J\f F[a, a^) is very useful in practice. It allows 
one to evaluate easily the expectation values of F{a, a''") in such important classes of 
states as the vacuum and coherent states, the latter defined for a complex number z 

n\n), {n\n') = Sn,n')j ^^^1 satisfying 



as \z) 



exp(-^)En=0 7^l'^), (with a^a\n) 
a\z) — z\z). 

Here we will apply the results of Section 2 and obtain J\f 
arbitrary functions q{x) and v{x). First note that 

'- fe=i 
with appropriately defined fn,k and Consequently one obtains 

gA[.(at)a+.(at)] ^ 1 + g [/.^(^t) + g /„,(at)a'=) , 

n=l \ k=l J 

whose matrix elements between the coherent states \z) and \z') are 

(^/|eA[,(at)a+.(at)]|^^ = {A^) (l + E^ ^ + E^ /nfe(^'*)^')) ■ (21) 

The next objective is to evaluate the sum on the r.h.s. of Eq.(21). We do so using the 
representation — > x and « — ^ and rewrite Eq.(20) as 



for 
(19) 

(20) 



= 1 + g a; [u^^^) + g !nk{x)^ 



n=l 



k=\ 



Next, by acting with this operator identity on e^^ we have 



gA[g(a:)£+«;(x)]gya; 



1+ E 



n=l 



hn{x) + E fnk{x)y^ 
k=l 



(22) 



(23) 



Eq.(l) allows one to rewrite the l.h.s. of the above equation (with T(A, x) and g{\^x) 
solutions of Eqs.(2) and (3)) as 

eM^)±+Mey^ ^ g{\ x)e^^(^'^), (24) 
which, if put into Eq.(23), implies 

oo yn 



1+ E 



n=l 



hn{x) + E fnk{x)y^ 
k=l 



(25) 



Using the last equality on the r.h.s of (21) we obtain 

(^/|eA[5(at)a+.(at)] _ ^_^,|^^ ^^^^ ^/*^g[T( (26) 

We now apply the crucial property of the coherent state representation (see [15] and 
[16]). This is that if for an arbitrary operator F{a, a^) we have 

{z'\F{a,a^)\z) = {z'\z) G{z'*,z) 
then the normally ordered form of F(a, a^) is given by 

Af[F{a,a^)] =: G{a^,a) : . 
Eqs.(26) and (28) then provide the central result 



(27) 
(28) 



^A[g(at)a+«(at)] 



=:^(A,at)e[^(^'"')-"1": , 



(29) 



being an operator identity in which the functions g, T are found through Eqs.(2) and 
(3). 
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4. Normal ordering and ShefFer-type polynomials 

We now investigate in more detail the properties of Eq.(26). In particular we shall 
identify its inherent characteristic polynomial structures. Using the representation of 
coherent states in terms of the Glauber displacement operator D{z) 

\z) = D{z)\0) = exp {za^ - z*a)\0) = e'l^l'/'e^'^' |0), (30) 

where |0) is the Fock vacuum, we rewrite the l.h.s of Eq.(26) as 

(^/|eA[9(at)a+.(at)]|^^ ^ 

^ g-l/2(k'P+|^|2)^Q|g^'.ag^atg_^„tgA[</(at)„+^(at)]^^„t|Q^_ ("^1) 

Next we apply the exponential mapping formula Eq.(5) in order to get 

g-zatgA[g(at)a+i;(ot)]g2at ^ ^\[q{a'*){a+z)+v{a'*)] 



,za' „z a„z z 



(32) 



which, if put into Eq. (31), give (with {z\z) = exp{z'*z — \\z'\^ — \\z\'^)) 

If we expand the exponential on the r.h.s. of Eq.(33) as a Taylor series in A, then the 
expansion coefficients for any fixed z' are polynomials in z. 

Properties of such polynomials may be clarified when one recalls the definition of 
the family of Sheffer-type polynomials Sn{z) [9]. The latter are defined through the 
exponential generating function (egf ) as 

oo \ Jl 

l + E5,(^)-=>l(A)e^^W (34) 

where functions A{\) and B{\) satisfy: A(0) = 1 and S(0) = 0, 5'(0) ^ 0. Many 
properties of these polynomials may be elucidated by means of the so called umbral 
calculus [9], [17], [18] which provides us with numerous interesting applications. 

Returning to normal order, recall that the coherent state expectation value of 
Eq.(29) is given by Eq.(26). When one fixes z' and takes A and z as indeterminates, 
then the r.h.s. of Eq.(26) may be read off as an egf of Sheffer-type polynomials defined 
by Eq.(34). The correspondence is given by 

AiX)^giX,z'*), (35) 
B{X) = [r(A, z'*) - z'*] . (36) 

This allows us to make the statement that the coherent state expectation value 

of the operator exp X{q{a^)a + v{a^)) for any fixed z' yields (up to the overlapping 

factor {z'\z)) the egf of a certain sequence of Sheffer-type polynomials in the variable z 
given by Eqs.(35) and (36). The above construction establishes the connection between 
the coherent state representation of the operator exp X{q{a^)a + v{a^)) and a family 
of Sheffer-type polynomials S^'^{z) related to q and v through 



\ n=l J 
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where explicitly (again for z' fixed): 

5^9.")(^) = {z'\z)-\z'\[q{a^)a + v(a^) 



\z) 



(38) 



e5l^'l'(z'| [q(a^)(a + z)+v(a'')Y\0). 

We observe that Eq.(38) is an extension of the seminal formula of Katriel [2], [6] where 
v{x) — and q{x) — x. The Sheffer-type polynomials are in this case exponential 
polynomials [9] expressible through the Stirling numbers of the second kind. 

Having established relations leading from the normal ordering problem to Sheffer- 
type polynomials we may consider the reverse approach. Indeed, it turns out that for any 
Sheffer-type sequence generated by A{X) and -B(A) one can find functions q{x) and v{x) 
such that the coherent state expectation value {z'\ exp \{q{a'')a + v{a'')) \z) results in 
a corresponding egf of Eq.(34) in indeterminates z and A (up to the overlapping factor 
{z'\z) and z' fixed). Appropriate formulas can be derived from Eqs.(35) and (36) by 
substitution into Eqs.(2) and (3): 

q{x) = B'{B-\x- z'*)), (39) 

, , A'(B-^x-z'*)) 
v(x) = — -i ^. (40 

One can check that this choice of q{x) and f (x), if inserted into Eqs. (2) and (3), results 
in 

T{\,x) = B{\ + B-^{x-z'*)) + z'\ (41) 

^^^'''^ - AiB-^ix-z'*)) ' ^^^^ 
which reproduce 

^^'|eA[g(at)„+.(at)]|^^ _ (z'|z)A(A)e^^(^). (43) 

The result summarized in Eqs. (35) and (36) and in their 'dual' forms Eqs. (39)- (42), 
provide us with a considerable flexibility in conceiving and analyzing a large number of 
examples. 



5. Combinatorial structures 



In this section we will work out examples illustrating how the egf (= S^o^(^)^) of 
certain combinatorial sequences a{n), n = 0, 1, 2, . . . appear naturally in the context of 
boson normal ordering. To this end we shall assume specific forms of q{x) and v{x) thus 
specifying the operator that we exponentiate. We then give solutions to Eqs. (2) and 
(3) and subsequently through Eqs. (35) and (36) we shall write the egf of combinatorial 
sequences whose interpretation will be given. 

a) Choose q{x) = x^, r > 1 (integer), v{x) = (which implies g{\,x) = 1). Then 
(see Eq.(15)) r(A,x) = x [1 - A(r - l)x^-^]^. This gives 

at \ 1 

: exp I J l a : (44) 

■ - A(r-l)(at)'-i)— J \ 



N 
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as the normally ordered form. We now take z — 1 in Eqs.(2) and (3) and from Eq.(26) 
one has 



exp 



(l|z)-Vl|e^('^')'''^|z 

\(1-A(r-1))- 
which for 2; = 1 generates the following sequences: 

r = 2 : a{n) = 1, 1, 3, 13, 73, 501, 451, ... 
r = 3 : a{n) = 1, 1, 4, 25, 211, 2236, 28471, . 

These sequences arc enumerating r-ary forests [3], [19], [20], [21]. 

b) For q{x) = xln(ea;) and v{x) = (implying g{X,x) = 1) which give T(A,x) 
e^^'^x^^ , this corresponds to 



etc. 



(45) 



(46) 
(47) 



^Xa^ \n{ea^)a 



: exp 



(48) 



(49) 



whose coherent state matrix element with z' = 1 is equal to 

(l|^)-^(l|e^"^^°(""^)"|^) = exp [z (e"^-^ - l)] , 

which for 2; = 1 generates a{n) — 1, 1, 3, 12, 60, 385, 2471, ... corresponding to partitions 
of partitions [19], [20], [21]. 

The following two examples will refer to the reverse procedure, see Eqs.(39)-(42). 

We choose first a Shcffcr-typc cgf and deduce q{x) and v{x) associated with it. 

c) A{X) = B{X) = A, sec Eq.(34). This egf for z — 1 counts the number of 
arrangements a{n) = n\ X]fc=o = li 2, 5, 65, 326, 1957, ... of the set of n elements [22]. 
The solutions of Eqs.(39) and (40) are: q{x) = 1 and v{x) 
corresponds to 



. In terms of bosons it 

2—x 



M 



A a+- 



2 - at - A 



2 - at - A 



(50) 



d) For A{X) = 1 and B{X) = 1 — \/l ^ 2A one gets the egf of the Bessel polynomials 
[23]. For z = 1 they enumerate special paths on a lattice [24]. The corresponding 
sequence is a(n) = 1,1,7,37,266,2431,... . The solutions of Eqs.(39) and (40) are: 
q{x) 



and v{x) = 0. It corresponds to 



A— 

e 2- 



^l--^(2-at)-2Aja 



(51) 



in the boson formalism. 

These examples show that any combinatorial structure which can be described by 
a Sheffer-type egf can be cast in boson language. This gives rise to a large number of 
formulas of the type Eqs.(50) and (51) which are important for physical applications. 
Many other examples can be worked out in detail and will be listed in a subsequent 
publication. 
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6. Conclusions 

The formula of Eq.(29) is the key result of the present investigation. It allows one 
to obtain normal ordering for a vast range of problems involving one annihilation 
operator (or one creation operator through hermitian conjugation). We have exploited 
the connection between the Eq.(29) and the egf of Sheffer-type polynomials and have 
given examples of combinatorial structures described by particular choices of q{x) and 
v{x). 

We thank G.A. Goldin for important discussions. 

PB wishes to thank the Polish Ministry of Scientific Research and Information 
Technology for support under Grant no: 1P03B 051 26. 
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